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Abstract 

We consider a cosmology in which a spherically symmetric large scale inho- 
mogeneous enhancement or a void are described by an inhomogeneous metric 
and Einstein's gravitational equations. For a flat matter dominated universe 
the inhomogeneous equations lead to luminosity distance and Hubble con- 
stant formulas that depend on the location of the observer. For a general 
inhomogeneous solution, it is possible for the deceleration parameter to differ 
significantly from the FLRW result. The deceleration parameter go can be 
interpreted as q$ > (qo = 1/2 for a flat matter dominated universe) in a 
FLRW universe and be qo < as inferred from the inhomogeneous enhance- 
ment that is embedded in a FLRW universe. A spatial volume averaging of 
local regions in the backward light cone has to be performed for the inhomo- 
geneous solution at late times to decide whether the decelerating parameter 
q can be negative for a positive energy condition. The CMB temperature 
fluctuations across the sky can be unevenly distributed in the northern and 
southern hemispheres in the inhomogeneous matter dominated solution, in 
agreement with the analysis of the WMAP power spectrum data by several 
authors. The model can possibly explain the anomalous alignment of the 
quadrupole and octopole moments observed in the WMAP data. 
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1 Introduction 

The problem of how to explain the accelerating expanding universe given the SNe 
la supernova data QjJ 0] in combination with the WMAP data has led to a host of 
solutions, ranging from modified gravity theories to a quintessence field and a cos- 
mological constant [3 H] . In the following we shall take the position that the homo- 
geneous and isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology is 



1 



the most symmetric model which describes many of the basic features of the data. 
However, important physical features may not be explained by this first approxima- 
tion and a more complete description is given by an inhomogeneous cosmology 
We shall investigate a universe in which a spherically symmetric inhomogeneous 
enhancement or void is embedded in an asymptotic FLRW universe. We find that 
the inhomogeneities permitted by observation can lead to a reinterpretation of the 
luminosity distance of a cosmological source in terms of its red shift z. The time 
evolution and the expansion rate of the inhomogeneous universe can lead to intrinsic 
effects such as cosmic variance at large angles, long-wavelength perturbations not 
described by a FLRW homogeneous and isotropic universe and late-time nonlinear 
effects. It is possible that the observed acceleration of the universe is caused by 
late-time inhomogeneity associated with large scale structure [0]. This is important, 
for it is difficult to explain theoretically the postulated dark energy that causes the 
acceleration of the universe. Tomita has investigated the interpretation of the lu- 
minosity distance and red shift in a local void [HE]. However, the determination 
of the sign of the deceleration parameter q depends on performing a spatial volume 
averaging of physical scalar quantities, such as the time evolution of the expansion 
parameter 9 in the Raychoudhuri equation [Oj. 

Recently, Barausse, Matarrese and Riotto, [HEI], Kolb et al. |IX| and Notari J2| 
have perturbed a homogeneous and isotropic FLRW universe to second order, and 
shown that the luminosity distance-red shift relation implies a non-vanishing cosmic 
variance of the deceleration parameter implying an intrinsic theoretical uncertainty 
in its determination. They uncover an infrared long-wavelength effect that depends 
on the fluctuation spectrum n s . For n s — 1 <C 0, a unit variance is obtained in second 
order perturbation that generates a significant cosmic variance in the deceleration 
parameter go for cosmological perturbations with the largest wavelengths. These 
results have been criticized by various authors [T3J EJ E3 EH H3 ITHIIT^] . 

Several of these investigations rely on using specific gauge conditions to carry 
out the perturbation calculations. Unless all the calculations of physical quantities 
are performed within a fully gauge invariant formalism, the choice of specific gauges 
such as the synchronous gauge or the Poisson (Newtonian) gauge can produce gauge 
dependent results. Moreover, some of the calculations assume an FLRW background 
spacetime. Rasanen [20] has considered the backreaction on a dust dominated 
Lemaitre-Tolman-Bondi (LTB) background spacetime 122 1 E ffi 12 ^ l2l ^ I27 H 15]. 
Perturbation calculations fail in the late-time nonlinear regime of large scale galaxy 
and void structure for red shifts z < 25. 

We shall not assume that the universe is a global FLRW model. Instead, we 
use an exact solution for the inhomogeneous enhancement, based on field equations 
for an inhomogeneous spherically symmetric distribution of matter, embedded in 
an asymptotic FLRW universe, specializing to the matter dominated LTB model. 
The assumption of spherical symmetry results in only one inhomogeneous spatial 
dimension. 

There have been recent reports of statistically significant anomalies in the cosmic 
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microwave background (CMB) compared to the standard big bang model [23 
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WMAP data appear to reveal anomalies at the largest angular scales > 60°. The 
angular two-point correlation function is suppressed at scales larger than 60°. When 
related to Fourier space, the vanishing of the two-point correlation function at large 
scales reveals an alignment of the quadrupole and octopole moments. This could 
be due to cosmic variance. Oliveira-Costa et al., [4T] studied the quadrupole and 
octopole data and showed that both are planar and aligned with all maxima and 
minima falling on a great circle on the sky. Similar results have been reported by 
Schwarz et al. [44] and Weeks [43]. It is not expected that these anomalies in the 
WMAP data are due to instrumental failure. The CMB temperature fluctuations 
have been found to be unevenly distributed between the southern and northern 
hemispheres with a statistical significance ~ 2 — 3a. Such effects could be produced 
by foreground contamination. If, however, they are true cosmic effects, then they 
would not agree with the standard predictions of inflation theory and a significant 
fine-tuning of inflationary potentials and parameters would be required to explain 
the phenomenon |49j . 

Schwarz et al., have applied a method assigning directions to the £-th multipole 
using multipole vectors [44] . They found that the method reveals a high statistical 
significance (99.9%) CL that the observed quadrupole and octupole are not in agree- 
ment with a Gaussian random, statistically isotropic sky as predicted generically by 
inflation. Surprisingly, they discovered a strong correlation with the orientation of 
the ecliptic plane (solar system and its motion) but no significant correlation with 
the Milky Way. The motion of the solar system is related to the measured CMB 
dipole. This effect could be due to a significant systematic error in the WMAP data 
or that the largest scales of the CMB sky are dominated by a local foreground. It 
could also be due to cosmic variance and suggest that we are in a special place in the 
universe at the present time. The cosmic variance for Gaussian random variables is 
given by 



It is a serious limitation for low multipoles that cannot be avoided. If, on the 
other hand, it is due to a true cosmic effect, or an absorbing or emitting source of 
microwave radiation in the solar system, then this will have significant consequences 
for interpreting cosmological data. 

Normally, the known Doppler shifts due to the motion of Earth round the Sun 
(30kms _1 ), the motion of the solar system in the Galaxy (~ 250 km s -1 ) and the 
motion of the Galaxy in the Local group of about 600 km s -1 in the direction I ~ 
260 °, b ~ 30° are removed from the data. Any residual velocities observed would not 
be expected to produce an alignment of the axes of lower multipoles. An alternative 
possibility is the existence of inhomogeneities on a scale approaching the scale of 
the visible universe itself. The CMB data show that such inhomogeneities must 
have a density contrast less than 10~ 4 over scales of order c/H . The deviation 
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of the observer from Hubble flow induced by the perturbations give a variation in 
temperature of amplitude (SB EE] : 
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where a s is the scale of the perturbation and q is the deceleration parameter for 
t = to and z = 0. The gravitational potential effects produce a fluctuation of order 
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where $ is the gravitational potential. The angular scale of the perturbation is of 
order a s H /c and at the surface of last scattering it is 



(Sp/p) r = (1 + z r )[a sr (l + z r )H r {l + z r )~^ 2 \' 



(4) 



where the subscript r denotes the quantity at the surface of last scattering (recom- 
bination). 

We have approximately 

— vdv, (5) 

for a perturbation of excess mass 5M with an induced peculiar velocity 5v in the 
Hubble flow v = HqR, where R is the distance from the center of the density 
perturbation. This corresponds to a density contrast 
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We shall investigate the interpretation of the red shift and luminosity distance 
relations in the LTB model. We will also determine the possible lack of complete 
smoothness of the temperature fluctuations and correlation functions across the sky 
predicted by the LTB model. The solution can be shown to have a high degree 
of isotropy and homogeneity at large scales, but there is a generic component of 
inhomogeneity and anisotropy that is direction dependent and correlated with the 
location of the origin in spherical polar coordinates. 



2 The Lemaitre-Tolman-Bondi Model 

The CMB temperature effects were first discussed by Sachs and Wolfe [20] and Kan- 
towski considering small density contrasts over large scales for an Einstein-de 
Sitter universe and that the observer is located in a perturbation, or the pertur- 
bations were so distant that the radiation passing through them occurred at last 
scattering. Other situations were discussed by Rees and Sciama [HI] and also by 
Dyer [SH]. Raine and Thomas [33] analyzed the situation in which the observer is 
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situated near the edge of a large-scale small amplitude density enhancement in an 
open universe. We shall investigate in the following, a small amplitude and large 
scale spherically symmetric density distribution with the observer in a spatially flat 
universe. 

The inhomogeneous line element in comoving coordinates can be written as (see 
Appendix A): 

ds 2 = dt 2 - R' 2 (t, r)f- 2 {r)dr 2 - R 2 (t, r)d{l 2 , (7) 

where we choose units G = c = 1, dVL 2 = d6 2 + sin 2 9d(p 2 and f(r) is an arbitrary 
function of r only. For the matter dominated LTB model with zero pressure p = 
and zero cosmo logical constant A = 0, the Einstein field equations demand that 
R(t, r) satisfies 

2RR 2 + 2R(1 - f) = F(r), (8) 

with F being an arbitrary function of class C 2 , R = dR/dt, and R' = dR/dr. 
There exist three possible solutions depending on whether f 2 < 1,= 1, > 1 and 
they correspond to elliptic (closed), parabolic (flat), and hyperbolic (open) cases, 
respectively. 

The proper density of matter can be expressed as 

F' 

P= IQ-kR'R 2 - (9) 
By using (JHJ), we can solve Q to obtain 

n-i s £-i = ^(Iz/!_ 2 /A (io) 

p c 3H 2 S \ R 2 RRT K ' 

where 

H 2 s = l(Hl + 2H ± H r ). (11) 

Here, we have defined two Hubble parameters H r (t,r) and H±(t,r) for the local 
expansion of a spherically symmetric density perturbation, corresponding to the 
radial and perpendicular directions of expansion, respectively. We have 

There exist three possibilities for the curvature of spacetime: 1) f 2 > 1 open 
(O - 1 < 0), 2) f 2 = 1 fiat (fi - 1 = 0), f 2 < 1 closed (fi - 1 > 0). We define a 
critical density in analogy with the FLRW model for f 2 — 1: 

R 2 R R' , 
8 ^=W + 2 RK- (13) 

This corresponds to the critical density for flat spatial sections t = const. We can 
obtain the analogy of the Friedmann equation for a spatially flat density perturba- 
tion: 

Hi, A (14) 
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The total mass of matter within comoving radius r is 

Mir) = - r drf^F 1 = 4tt f drpf^R ' R 2 . (15) 
4 jo Jo 

Since the WMAP data [53] shows that the universe is spatially flat to within a 
few percent, we shall consider the globally flat case f 2 — 1. The metric reduces to 

ds 2 = dt 2 - R' 2 dr 2 - R 2 dVl 2 , (16) 

where 

R(t,r) = r[t + f3(r)] 2/3 , (17) 
and f3{r) is an arbitrary function of r of class C 2 [24 . The metric (fTo^) becomes 

ds 2 = dt 2 - (t + (3fl\Y 2 dr 2 + rW), (18) 

where 

2r8' 

Y ^ + WWY (19) 

and 

p = e^tTWr - m 

The arbitrary function /3(r) can be specified in terms of a density on some space- 
like hypersurface t = t . The metric and density are singular on the two hypersur- 
faces t + (5 — and Y = 0, namely, ti = —(3 and t 2 = — (3 — 2rj3'/3, respectively. 
The model is only valid for t > S(r) = Max [t\ (r ), t 2 (?")], and the hypersurface 
t(r) = £(r) defines the big-bang. However, our pressureless model requires that 
the surface t(r) = E(r) describes the surface on which the universe becomes matter 
dominated (in the FLRW model this occurs at z ~ 10 4 ). We observe that even in 
the spatially flat LTB model, different parts of the universe can enter the matter 
dominated era at different times. 

For (3 = and in the limit t — > oo we obtain the Einstein-de Sitter universe 

ds 2 = dt 2 -a 2 (t)(dr 2 + r 2 dn 2 ), (21) 

where a(t) = t 2//3 . Thus, for (3 = we obtain the FLRW model. Moreover, the 
expanding flat LTB model necessarily evolves to the homogeneous and isotropic 
FLRW model for a non-vanishing density, whatever the initial conditions. 



3 Paths of Light Rays in the Inhomogeneous Cos- 
mology 

The luminosity distance between an observer at the origin of coordinate system to, 
and the source at (t e , r e ,9 e ,(p e ) is 

/ £ \l/2 

d L = =R(t e ,r e )[l + z(t e ,r e )] 2 , (22) 
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where C is the absolute luminosity of the source, T is the measured flux and z(t e , r e ) 
is the red shift (blue shift) for a light ray emitted at (t e ,r e ) and observed at (to,0). 

For 6 = (ft = constant and f 2 (r) = 1, we have for a light ray travelling inwards 
towards the center and X(r,t) = R'(r,t): 

^ = -R'(r,t), (23) 

where the sign is determined by the fact that the light ray travels along the past 
light cone. Consider two rays emitted by a source with a small time separation r. 
The equation of the first ray is 

t = T(r), (24) 

while for the second ray, we have 

t = T(r)+T(r). (25) 
The equation of a ray and the rate of change of r(r) along the path is 

dT{r) 



dr 
dr(r) 



= -R'[T(r),r], (26) 
-r(r)R'[T(r),rl (27) 



dr 
where 

R'[T(r),r}=R'\ r , T{r) . (28) 
Choosing r(r e ) to be the period of a spectral line at r e , we get 

4^ = 4Sr = 1 + ^)' z = for r e = °- ( 29 ) 

r(r e ) i/(0) 

The red shift considered as a function of r along the light cone is determined by 

dz 

^- = (l + z )R'[T(r),r\. (30) 

For a light ray travelling from (ti,ri) to (to, 0) the shift Z\ is 

ln(l + Zl ) = - r drR'[T(r),r}. (31) 
Jo 

We can compare this to the result obtained from a FLRW universe model by re- 
labelling the radial coordinate f = R[T(r),r] and choosing the solution f 2 (r) = 1, 
which corresponds for p = A = to a flat, inhomogeneous matter dominated solu- 
tion. Moreover, we choose 

ai (r) = R[T(r),r], a 2 (r) = R[T(r),r], ... (32) 
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as functions of r only. By expanding (j31j) . we obtain 

ln(l + «i)=/ dr— -± / dr^ = -ln(l - ai ) - / dr ., V \ . (33) 

Jo 1 — Oi jo 1 — a\ Jo (1 — ai) 

Two terms contribute to the cosmo logical red shift: the contribution due to the 
expansion of the universe, and the shift due to the difference between the potential 
energy per unit mass at the source and at the observer. For the FLRW case, M'(r) = 
and there is no gravitational shift. The integral in Eq. (j33)l for small r\ can be 
neglected, and expanding the logarithms on both sides we find for small r e or small 
to — t e : 

z(t e , r e ) = H ± (t e , r e )d L (t e , r e ). (34) 

The difference between this result and the one obtained from FLRW is that it is 
local and depends on the angular Hubble parameter H± = R/R rather than on the 
FLRW Hubble parameter Hflrw = d/a. 

A formula for the luminosity distance in an exact FLRW universe is given by 



(d L ) FLRW (z) = c(l + z) ~Tj~7 \ • (35) 

JO n[u) 



This can be reexpressed as 

<•* dz 



(d L ) FLRW (z) = c(l + z)\l - Qol^H^S |1 - Q \ 1/2 H l Z 

L Jo 



Hz) 



(36) 



where S(x) = sin(a;) for Qq > 1, sinh(x) for Q < 1, and x for Q Q = 1. For a flat 
universe 



(d L ) FL Rw(z) =c{l + z)H 1 [ Z duexp - Hi + q(v(r))]d\n(l + v) 

jo L Jo 



(37) 



where q denotes the deceleration parameter. This result only depends on the as- 
sumption of a FLRW universe and does not depend directly on solutions of Einstein's 
gravitational equations. In our inhomogeneous cosmology, we must generalize the 
formula for d^. 



{d L ) LTB (z(r)) = c(l + z(r))Hll [** duexp\- Hi + q(v, r)]dln(l + v 

Jo Jo 



(38) 



Thus, we interpret the angular Hubble parameter H±(t,r) = R(t,r)/ R(t,r) as the 
Hubble parameter that replaces the FLRW expression H(t) = a(t)/a(t). 
Let us expand R(r, t) in a Taylor series 

R(r, t) = R{r, t Q - (t - t)} = R(r, t )[l - (t - t)^44 + J(*o - t? ^ ^ 1 



R(r,t ) 2 K 1 R(r,t ) 
R(r,t )[l - (t - t)Ho± - i(t - t) 2 g (r,to)^o\ " ■■■), (39) 
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where to denotes the present epoch and H Q ± = R(r,to)/ R(r,t ). Moreover, we have 



5o(rto) = _flM^o) } 

R 2 (r, t ) 



By substituting for R from Eq. (j81|) (Appendix A), we obtain 

1 47rp (r, t ) _ A _ 1 H 0r {r,to) 
••5 ' ^/At. (r./ ) 3tf \Mo) 3U oi (M„)' 



where #or(Mo) = #(r, t )/R'{r, t ). 

If we set H ±(r,to) = Ho r (r,t ) = H where H = a(to)/a(t ), then we obtain 
the spatially flat FLRW expression for the deceleration parameter: 



^Po A 1 
q ° = MlJ~ 3H$ = 2 QoM ~ QoA - (i2) 



We see from (|41|) that different observers located in different causally disconnected 
parts of the sky will observe different values for the deceleration parameter q Q , de- 
pending upon their location and distance from the center of the spherically sym- 
metric inhomogeneous enhancement. 

We must interpret the Hubble parameter and red shift measurements in our 
model given a void on the scale of ~ 100 Mpc or a large scale inhomogeneous 
enhancement of the scale of the Hubble horizon ~ c/H. Observers at spatially 
separated locations in the universe would interpret measurements of z, H and di 
differently from one another. We interpret the deceleration parameter in the FLRW 
model as 

% = ~ n(t (43) 

a{t )li 

and from the SNe la supernovae observations arrive at a value go < corresponding 
to an accelerating universe. However, an observer in the inhomogeneous model 
universe in a different, causally disconnected location in the universe could observe 
a different go an d luminosity distance d^. Thus, averaging all the local observations 
of di and evaluations of go can l ea d to a cosmic variance on the determination of 
the deceleration parameter. A significant generic, theoretical uncertainty (cosmic 
variance) is generated in the determination of the deceleration parameter go. 

Luminosity distances are determined by supernovae measurements, which cannot 
directly determine the instantaneous expansion rate or deceleration rate. Thus, to 
determine the expansion history, assumptions must be made about the evolution 
of H(z) and q(z). We find that our inhomogeneous, spherically symmetric model 
of large scale inhomogeneous enhancements, reaching out to the Hubble horizon 
scale c/H, can lead to significant differences in the local determinations of H and 
di and an observer independent inference of the rate of expansion and deceleration 
of the universe has a generic, built-in theoretical uncertainty. If the universe were 
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locally LTB for a matter dominated universe rather than FLRW, then the Hubble 
parameter based on the observed LTB values of z determined through the relation 
(|38|) would be position and red shift dependent. The dependence of H versus z 
and z versus di were given in [SJ 122] for different evolving ages of the universe. It 
was found that H(z) strongly depends on the choice of model for the FLRW case, 
and the choice of the deceleration parameter qo in pTj) and (|38j) through which we 
interpret the observations. The calculations in [01 ESj were based on the equation 
describing the ratio of the local density to the density of the FLRW universe taken 
along the past light cone for a void: 

V ° ld[ 1 )l ~ p FLR w[T(r)] 8R>[T(r),r]R*[T(r),rY W 
where j3 is a constant. 



4 CMB Fluctuations and Inhomogeneity 

The main assumptions of the standard big-bang model are that the universe is 
spatially flat and dominated in the matter era by cold dark matter (CDM) and 
dark energy. The universe enters the matter dominated era with scale invariant, 
Gaussian, adiabatic initial fluctuations generated in the early universe by a period 
of inflation. Through a hierarchical gravitational instability process these uniform 
fluctuations started growing after decoupling (surface of last scattering) to form the 
large scale structure in the form of galaxies and clusters of galaxies that we see 
today. 

In spatially flat models, the density contrast is define by 

i(x) = 6 -M = PW-P, (45) 
P P 

where p denotes the mean density. Its Fourier transform is 

6k=yJ d 3 x6(x) exp(ik • x). (46) 

The power spectrum \5k\ 2 corresponds to statistical random gaussian fluctuations 
with the rms density fluctuations 

^ = (5(x)5(x)>V2. (47) 
P 

For an isotropic power spectrum (i.e. depending on k = |k| instead of k): 

Sp 1 f°° k 3 \S k \ 2 dk , , 

/ 1 ■' (48) 



p V Jo 2tt 2 k 
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The autocorrelation function is defined by 



e(r) = <<y(x + r)<5(x)>, 



(49) 



which is the Fourier transform of the power spectrum \5k\ 2 and £(0) = (5p/p) 2 . 

Since our LTB model is confined to the matter dominated era, it is reasonable to 
assume that to a first approximation in our inhomogeneous flat model the density 
perturbations grow like 5 + = (Sp/p) oc t 2 / 3 independently of the scales. Moreover, 
we assume that the fluctuations entering the matter dominated era on the surface 
t(r) = S(r) have the inflationary scale invariant behavior of the CDM FLRW model. 
In the FLRW model, we have for the growing mode [2*K] : 



$FLRw(t) — ^FLRwiteg) 



ft FLRW \ 2 / 3 



V t 



eq 



(50) 



where t eq and tpiRW denote the time of decoupling and the time from the initial 
FLRW singularity to a given coordinate value of time t, respectively. In the FLRW 
model these times are the same everywhere in space. 

In the LTB model for a spatially flat universe, we have 



'LTB 



{t,r) 



'FLRW 



tLTB{r) 



2/3 



(51) 



where now Iltb{t) = t — ts(r) is the time from the initial singularity t > ts(r) to a 
given coordinate time t. For some time coordinate value t: 



5 L TB{t,r) = S FLRW (t) 
We have made the simplifying assumption: 



tLTBjr) 
tFLRW 



2/3 



8FLRw(teq) = hTB(tT,(r),r) and t eq = tx(r) forallr. 



(52) 



(53) 



For some spacelike hypersurface t = const, the gravitational amplification of 
the primeval fluctuation perturbations depends on the spatial position. The larger 
^LTs( r ) for a given r, the more developed structure we expect to see. Also, there 
can be a different observable distribution of fluctuations in different parts of the sky 
depending on t LTB {r) and the location of the observer. 

We can now write the two-point correlation function as 



fz,TB(*o,r) = (hTB(to,r)S LTB (t ,0)). 
The correction to the FLRW correlation function ^flrw is given by 

"Wr)l 2 / 3 



tFLRW 



(54) 



(55) 
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The CMB temperature anisotropies generated by scalar perturbations of the 
FLRW spacetime metric can be written as 



ds 2 = a 2 (r])[(l + 2§)drf - (1 - 2§)dx 2 



(56) 



where a{rj) is the FLRW scale factor and 77 is the conformal time. At the last 
scattering hypersurface, the anisotropy contributes $/3 sources. For a non-zero 
d$ / drj an additional integral contribution arises, for a scale free primordial spectrum 
|$ fc | 2 = \A 2 /k 3 \ the CMB angular power spectrum amplitude can be written 



where 



6T\* 



^(li±l T \A 2 K 2 



\£(£ + l) 



1 1 



(57) 



r°° rile rvo riti 

Kj = 2£(£ + 1) j o ™\j e (kVo) + 6 J d v ^j e (k( Vo - r?))] 2 , (58) 

and dfk/dr) denotes d^k/dr] normalized to A/k 3 ^ 2 . A flat CDM model has dQ/dr/ = 
and K 2 = 1, while a flat AC DM has K 2 > 1 with an amplification exceeding the 
Sachs-Wolfe plateau. 

For the late time universe, it was argued by Contaldi et al. jlH] that the sup- 
pression of the large scale (low £ moments) power spectrum could be explained by 
fine-tuning the coefficients K 2 to be small with certain solutions of late time $. On 
the other hand, the amplitude of metric perturbations at horizon crossing in infla- 
tionary models could fine-tune fc 3//2 $fc to give a small contribution at large angular 
scales. This would involve fine-tuning the shapes of inflaton potentials in an ad hoc 
way. 

In our LTB model, the growth of the angular power spectrum will have the form 



ST 



T ))F B {t,r) = ({^y t L ™{t) 



tLTBir)] 2 / 3 



ty 



(59) 



We see that there will be a difference in the observed CMB temperature fluctuations 
in different parts of the universe, depending on the value of r and the location of 
the observer. 



5 Isotropy of the Cosmic Microwave Background 

The standard description of cosmology is based on the FLRW model and the cos- 
mological principle - the statement that the universe is the same in all directions 
and locations. The assumption of isotropy necessarily implies that the statistical 
properties of the CMB should be the same in all directions on the sky. The FIRAS 
experiment on the COBE satellite j^Hl IHZj demonstrated that the mean temperature 
of the CMB is isotropic to a high precision. However, the sensitivity of the WMAP 
data [S3] can be used to test the isotropy of the angular fluctuations in the CMB. 
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Hansen, Banday and Gorski have studied the WMAP power spectrum ex- 

tracted from the CMB using patches at different directions on the sky. In the lower 
angular range £ ~ 2 — 40, they find preliminary indications for non-zero differences 
in the power spectrum in the northern and southern hemispheres oriented along the 
galactic colatitude and longitude (80°, 57°), close to the ecliptic pole. Eriksen et 
al. j2H] and Park [21] were the first to discover a difference between the northern 
and southern hemispheres in non-gaussianity tests. Similar results have been found 
by Coles et al. (33 3 , Cruz et al. [Slj, Vielva et al. [SSj, Komatsu, Spergel and Wan- 
delt 36J, Larson and Wandelt [HZ], Land and Magueijo Oliveira et al. [JT]. An 
analysis of the variation in cosmological parameters associated with the acoustical 
peak in the power spectrum was carried out by Donoghue and Donoghue, to see 
whether there exists a correlation between the height and location of the peak [3Zj ■ 
A measurable dipole effect in the CMB would indicate that there exists a spatial 
asymmetry. However, at present after the subtraction of the dipole contribution 
from the WMAP data, there is no indication of a statistically significant dipole 
contribution. 

The observed non-zero differences in the power spectrum in the northern and 
southern hemispheres can be ascribed in our LTB model as being due to a difference 
in tz,Ts( r ) i n the locations of the northern and southern hemispheres, which will 
lead to a correction factor described in Eqs.(jH2J) an d (foTlj) : 



where r s and r n denote the position locations of the southern and northern hemi- 
spheres, respectively. 

The claims about the detection of an asymmetric distribution of large scale power 
in the CMB anisotropy as measured by the WMAP satellite are found at large an- 
gular scales. For scales £ > 40 (corresponding to an angular scale ~ 3 — 5°), the 
distribution of the CMB fluctuations is consistent with the hypothesis of isotropy 
and homogeneity, although there have been indications of possible foreground con- 
tamination around the first peak £ ~ 220. However, for lower £ multipole values, a 
strong difference between the northern and southern hemispheres (for galactic and 
ecliptic frames of reference) was found. The CMB power spectrum at large scales 
(£ = 2 — 40) is found to be significantly much lower in the northern hemisphere than 
in the southern hemisphere, which leads to a different estimate of the cosmological 
parameters in these hemispheres. 

When attention is focussed on three parameters to which the analysis of the 
anisotropy and non-Gaussianity is most sensitive, namely, the spectral index n s , the 
amplitude of fluctuations A and the optical depth r, the preferred value for the 
optical depth in the north is r < 0.08, whereas in the south it is r = 0.241q;q7 (68% 
confidence level). The latter result is inconsistent with r = at the 2a level. The 
WMAP collaboration estimate r = 0.17 could thus originate in structure associated 
with the southern hemisphere. By setting a prior on r, values of the spectral index 




(60) 
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n s are found that are inconsistent between the opposite hemispheres. 

Let us now assume that a mechanism is operative such that the amplitude (JHHJ) 
is suppressed for large angular scales in the sky for low multipoles. The observer is 
then able to see small dipole, quadrupole and octopole contributions generated by 
an inhomogeneous perturbation enhancement at scales c/H and the gravitational 
potential at the location of the enhancement. An observer at (ro,to) measures a 
temperature T = To. At the time of recombination z ~ 10 3 , the temperature is 
T r ~ 10 3 To. The red shift varies somewhat with angle 9 between the direction of a 
light ray emitted at (r ,t ) and absorbed at (r r ,t r ) and a light ray directed toward 
the center of the density perturbation. The apparent temperature of the CMB as 
observed in the direction 9 is given by 

™ = TTWr T "TTWY (61) 

where T av and z av are the temperature and red shift averaged over the whole sky 
and we have T av ~ T . 

The dipole, quadrupole and octopole moments, T>, Q and O are defined by 



V=^r r T b (9)Y l0 (9)sm9d9, (62) 

4v JO 

Q = J- rT b (9)Y 20 (9)sm9d9, (63) 
4v Jo 

= ^- T b (9)Y-w{9)sm9d9, (64) 

4v JO 

where e.g. Y w = ^j3/4n cos 9. An observer looking toward the center of the spheri- 
cally symmetric LTB density perturbation will see an axially symmetric distribution. 
The dipole V, quadrupole Q and the octopole O will be aligned in this plane, with 
the angle 9 defining the trajectory of a light ray arriving to the observer located 
at (r ,t ). The dipole, quadrupole and octopole contributions are produced by gra- 
dients of the gravitational potential associated with the large scale perturbation 
enhancement. The dipole contribution is small and undetected in the WMAP data, 
while the quadrupole and octopole contributions have been found to be planar and 
the planes are aligned to a statistically anomalous degree [JTJ 113 GS] • 

A possible explanation for the quadrupole and octopole anomaly is that the 
cosmic topology has the form of a toroidal universe with one small dimension of 
order one-half the horizontal scale, in the direction towards Virgo (small universe 
model). Another possibility is that the universe takes on the topological shape of a 
dodecahedron [60J. These possible explanations appear to have been ruled out by 
Cornish and collaborators |6*T] . 
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6 Conclusions 



We have used a cosmological model with the exact inhomogeneous LTB solution 
for a matter dominated, large scale inhomogeneous enhancement at Hubble horizon 
scale cj H with p = A = to demonstrate that observers located in different places 
in the universe can differ significantly in their determinations of the evolution of 
red shift and Hubble expansion rate. This will significantly influence the theoretical 
interpretation of the locally determined values of H(z) and d^z). Averaging over 
all the observers' results will yield a generic cosmic variance in the determination of 
the deceleration parameter. This in turn will lead to a variance in the conclusion 
as to whether the universe is undergoing an accelerating phase. An independent 
way to determine luminosity distances to supernova and other distant objects in the 
universe could settle the theoretical ambiguity. 

The strong difference between the distribution of CMB fluctuations in the north- 
ern and southern hemispheres discovered by analysis of the WMAP power spectrum 
data for large scales with multipole £ = 2 — 40 values, in the ecliptic and galactic 
frames of reference, can be explained as a cosmological effect by a large scale inhomo- 
geneous enhancement, described by our exact inhomogeneous solution of Einstein's 
field equations. The time evolution of the universe as measured from the surface of 
last scattering in the matter dominated LTB model can produce a correction factor 
(£ltb(t , )Ae) 2//3 for the fluctuations that can differ significantly between the northern 
and southern hemispheres. However, the uneven distribution of fluctuations could 
be due to systematic effects in the WMAP data or be due to foreground contamina- 
tion. Hopefully, the WMAP2 analysis of the CMB data will decide whether or not 
the non-Gaussian non-isotropic effect has a cosmological origin or is due to some 
other non-cosmological mechanism. 

The possible statistically significant alignment of the quadrupole and octopole 
moments in the WMAP data can be explained, in our large scale inhomogeneous 
model, by an off-center observer seeing an axisymmetric alignment of the quadrupole 
and octopole moments as the observer receives light rays from the center of the large 
scale inhomogeneous enhancement. 

7 Appendix A 

For the sake of notational clarity, we write the FLRW line element 

ds i = d1 ? _ a 2( t ) ( + r i dn A , (65) 
V 1 — kr 2 J 

where dfl 2 = d9 2 + sin 6 2 d<f) 2 and k has its usual interpretation. 

Let us now consider the more general, spherically symmetric inhomogeneous line 
element [2U 1221 123]: 

ds 2 = dt 2 - X 2 (r, t)dr 2 - R 2 (r, t)dVl 2 . (66) 
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The energy-momentum tensor T^ v takes the form 

T" v = {p+p^Uv-pS^, (67) 

where = dx^/ds and, in general, the density p = p(r,t) and the pressure p = 
p(r,t) depend on both r and t. We have for comoving coordinates u° = l,u l = 
(i = 1,2,3). Moreover, the total mass is 

M(r) = 4tt [ r T° ov ^ = 4vr f drp{r)XR 2 , (68) 
Jo Jo 

where g = Det(<7 Mi ,) = — X 2 R A sin 2 6. From this follows that 

M > = = AnpXR 2 . (69) 
dr 

We have assumed, in (|68|). that geodesies of particles do not intersect so that the 
total mass M = M(r) depends on r only [25] , 
The Christoffel symbols are 

r V = ^9 aa (d^g ua + d v g m - d a g^ u ). (70) 
The non-vanishing Christoffel symbols are 

pl _ X o _ • x _ X 2 _ 3 _ _R 

1 11 — "^r, J- 11— AA, 1 oi — J- 02 — J- 03 — 

r° 22 = -R-R, r 2 i2 = r 3 i 3 = — , r 1 22 = — r° 33 = rr sin 2 0, 

^33 = sin2 ^ r2 33 = - sin 9 cos 5, r 3 23 = cot B, (71) 

where X' = dX/dr and X = dX/dt. The Einstein gravitational equations are 

G^ + Ag^ = -87iT^, (72) 

where G^ u = R^ u — \g^ v R, R = g^R^ and A is the cosmological constant. We 
have 

G ° = 2^ + - i- ( 2^ + ^ - 2^) = SvrTo + A = 8vrp + A, (73) 

XR R 2 X 2 \ R R 2 XR) H ' V y 

Gi 1 - 2| + - j^ii = ^ + A = + A > ( 74 ) 

_ o X R XR 1 / i? X'R'\ „ o „ o ._ „ 
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Gi° = -2 - = SttT, = -SttXW = 0. (76) 

From Eq. (|76|) . we find that 

X(r,t) = -^-R'(r,t). (77) 
f( r ) 

For an isotropic pressure T^ 1 = T 2 = T 3 3 and the pressure p = pit) only depends 
on the time t. 

We now obtain the two equations 



i? , IR 2 1 1 „ x2N IR'R t Iff 4tt,_ , , , 1 
We obtain for f 2 = l the result 



H + 3W + 3 - ; > - 3 a r + 3 mi = "T (p + 3p) + 3 A ' (79) 



| +2 || = 8 . p + A , (80 ) 

* + - * + 3p) + iA. (81) 

R 3R 2 3R' R 3 KH y> 3 K ' 

By using the notation H± = R/R and H r = R'/R', we can write (JHUj) and ()81|) as 

Hi + 2H ± H r = 8np + A, (82) 
R 1 1 Att 1 

- + -Hi - -H r H ± = -—{p + 3p) + -A. (83) 

For H±(r,t) = H r (r,t) = H(t) = a/a and R(r,t) = a(t), we obtain the Friedmann 
equations of FLRW for a spatially flat universe: 

H 2 = ^ + 1 A. (84) 
3 3 V ; 

a 47T , . 1 . . 
- = - — (p + 3p) + -A, 85 
a 3 3 

where p = p(t) and p = pit). 
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